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Abstract. This paper is concerned with the cyclicity of the period annulus 
of quadratic Hamiltonian triangle vector field under quadratic perturbations. 
This problem has been studied by Iliev (J. Differential Equations 128(f996)), 
based on the displacement function obtained by Zoladek (J. Differential Equa- 
tions 109(1994)). Recently, P. Mardesic etc. (J. Dynamical and Control Sys- 
tems 17(2011)) studied unfoldings of the Hamiltonian triangle within quadratic 
vector fields. It turned out that the displacement function is not precise of the 
form given by Zoladek. Using the corrected form of the displacement function 
obtained by P. Mardesic etc, it is proved in this paper that the cyclicity of the 
period annulus under quadratic perturbations is equal to three. 



1. Introduction and the main result 

In this paper we study the bifurcations of limit cycles in a class of planar qua- 
dratic integrable systems under small quadratic perturbations. Taking a complex 
coordinate z = x + iy and using the terminology from [22] . the list of quadratic 
centers at z = looks as follows: 

z = —iz — z 2 + 2\z\ 2 + (b + ic)z 2 , Hamiltonian (Q^), 

z = —iz + az 2 + 2\z\ 2 + bz 2 , reversible (C}f), 

z = —iz + Az 2 + 2\z\ 2 + (6 + ic)z 2 , \b + ic\ = 2, codimension four (Q4), 

z = —iz + z 2 + (6 + ic)z 2 , generalized Lotka — Volterra (Q^ v ), 

z = —iz + z 2 , Hamiltonian triangle, 

where a, b, and c stand for arbitrary real constants. 

Recently, the authors of the paper [T7| studied unfoldings of quadratic Hamil- 
tonian triangle vector field within the quadratic vector fields. These unfolding 
have been studied by Zoladek [211 122] , Iliev [9j [10] and others. It turned out that 
the displacement function is not precise of the form given by Zoladek. 

Based on the displacement function obtained by Zoladek [22J, Iliev [9 J studied 
the cyclicity of the period annulus of quadratic Hamiltonian triangle vector field 
under quadratic perturbations. Due to the error in Zoladek's paper, the problem 
of the cyclicity of the period annulus of quadratic Hamiltonian triangle vector 
filed is still open. As usual, we use the notion of cyclicity for the total number of 
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limit cycles which can emerge from a configuration of trajectories (center, period 
annulus, a singular loop) under a perturbation. 

To state the main results of the paper [T7] , consider the quadratic system 



(1) X 



x = x[f3-(3x-((3 + l)y], 
y = y[— a + (a + l)x + ay), 



which has a first integral x a y l3 (l — x — y) with the integrating factor x a ~ l y l3 ~ l . 
System (CQ) has a center at p a p(a/(a + f3 + l),f3/(a + f3 + 1)). The vector field 
X is Darboux integrable. If a = /3, then X e Q%. If a = (3 = 1, then X is a 
Hamiltonian triangle system. 

According to [22J , the following system gives the general unfolding of the Hamil- 
tonian triangle system within quadratic systems: 



(2) X £ 



x = x[f3 — f3x — ((3 + l)y] + e x 2 + Eiy 2 , 
y = y[—a + (a + l)x + ay] + e y 2 + e 2 x 2 , 



where eoj^i an d £2 are small parameters. 

It is well known that the displacement function 5 is the difference between the 
first return Poincare map and the identity. We study the displacement function 
only in a compact domain within the basin of the center p a p away from the 
separatrices bounding the basin and the center itself. Set 

£ 3 = a — 1, £4 = P — 1, e = (eq, 61,62, 63, 64). 

The main result of the paper [T7] is the following: 

Theorem 1. [17] Let 

H £3 , E4 = x 1+£ *y 1+£ *(l-x-y). 
The displacement function S(h,e) of system ([2]) ; fore close to is of the form 
(3) 5{h, e) = /xi[Ji + •••]+ ix 2 [J 2 + ■■■} + /x 3 [J 3 + ■■■]+ fi 4 [J 4 + ■■•], 
where 

Ml = — £ 0) 

Hz = -(5^3 + e 2 e A )/2, 

fi 3 = (eiE 3 - e 2 e 4 )(e 3 - e A )/2, 



ji 4 = 



(£l£ 3 - E-iEij^Ex +£ 2 )/6, 
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and 



Ji — <£> y 2 dx — x 2 dy — —2 1 I (x + y)dx A dy, 

Jh=H 0}0 J Jh<H 0fi 



J 2 = / / (xy) x (x 3 + y 3 )dx A dy, 

J Jh<H 0fi 

J 3 = (xyY l [(x -y)(x + y) 2 \n(x/y)]dx A dy 

J Jh<H 0fi 

+ / / (xy)~ 1 (x + y)(x 2 + xy + y 2 )dx A dy, 

J Jh<H 0fi 

1 f x 3 y 3 

J ± = or f 7—z( dx + d y) 

on Jh=H . V x 

1 f f x 2 y 2 (3y 2 - Axy + 3x 2 ) , 
= oT / / 7 V2 dxAdy, 

3n J Jh<H 0)0 {x-yr 

and the three dots ■ ■ ■ denote analytic functions in (h,e), which tend to zero as 

As mentioned in [17] . the function J${K) does not vanish at the value of the 
center p a p, whereas Zoladek's does. This might seem surprising. An explanation 
is given in [T7] . 

In this paper we study the cyclicity of the period annulus of quadratic Hamil- 
tonian triangle vector field. The main result of the present paper is the following. 

Theorem 2. For small e, the maximum number of limit cycles in which 
emerge from the period annulus of system (JT]) is equal to three. The upper bound 
is sharp. 

Theorem [2] is proved by counting the number of zeros of the principal part J(h) 
of the displacement function S(e,h). To do it, we show that J'(h) has at most 
three zeros. Note that J(h) does not vanish identically at the value of the center 
Pa/3, whereas most of the principal part of the displacement function for other 
integrable systems do. Hence we need a careful analysis to get the exact upper 
bound of the number of zeros of J(h). 

Let us list here some results concerning the cyclicity of the period annulus for 
quadratic centers. The cyclicity of the period annulus for quadratic Hamiltonian 
Qf was completely solved by several authors, see [H [21 HI [HI [131 EE] and references 
therein. The generalized Lotka-Volterra was studied by Zoladek in [22J. 
Some results concerned with certain specific cases of Qf can be found in [3l [T5l[T6] . 
etc. The cyclicity of the period annulus for Q A was investigated in [51 [TS]. For 
more results on this problem, we refer the readers to the survey paper 



The rest of this paper is organized as follows. In Section 2 we give alternative 
expressions for the integrals J(h) and some preliminary results. Section 3 is 
devoted to prove the monotonicity and convexity of the ratio of two integrals 
w(h) = I' 2 {h) / I' Q (h) . The main result is proved in Section 4. 
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2. Preliminary results 

To use the results obtained in [9], we begin our investigation by a technical 
lemma giving alternative expression for the integrals Jk{h), k = 1, 2, 3, 4. 

Lemma 3. Let 

h(h) — f x l ydx, % — 1, 2, • • • , I*(h) = (b y(x — 1) \nxdx, 

where Th is the closed component of the level set 

(4) H(x, y) = x(y 2 - (x - 3) 2 ) = h, h e (-4, 0). 

Suppose that the oval Th has the positive (counterclockwise) orientation. Then 
the following relations between the integral 4 and 4 hold: 

2 

Ji = 27/0, 

4 = ^((/i + 18)J -12/ 2 ), 

4 = ^((/i- 12)4 + 244 -364), 

J4 = l^ ((19/i + 7 ° 2)/ °- 324/2) 

+ — ^—((-3888 - 324/i + 7h 2 )I' + 216(6 + h)I'). 
20736/i ; V ' 2J 

Proof. First of all we list several relations between Ik(h), which will be used in 
the proof of this lemma. It has been proved that the following two identities 
holds in [9]: 

(5) I (h) = h(h), (2k + 6)/ fc+ i = 6(2A; + 3)4 - 18*4-1 - (2* - 3)^4-2- 
One obtain from (j3j) that 

dH = 2xydy + (-9 + I2x - 3x 2 + y 2 )dx. 

Multiplying both sides of the above equation by x % yi~ 2 and integrating over 14 
yields 

= 2/ x i+l y j ~ l dy - 9/^-2 + 12J i +i )j _ 2 - 3/ i+2j _ 2 + 4, 



where 



On the other hand, 



Iij = Iij(h) = (b x l y 3 dx. 



.'•'•'// l dy = - — I„. 

J 
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Therefore, 

j — 2i — 2 

(6) Iij = 9Zjj_2 — 12/j + i > j_2 + 3/j+2J-2- 

It follows from fl4]) that 

(7) rf ( ln |^) = ^ lfa _?£(£^) d!/ , 

\ 3 — x + y J H H 

and 

(8) I'tih) = I x'%dx = <£ ^— dx. 



dh J Vh 2y 



Perform a change of the variables which transforms the Hamiltonian ifo,o into 
H, defined in (jlj). This is given by (x, y, iJo,o> h) — > (xi, yi, H, hi), where 

(9) x = \(3 - x x + yi), y=\{$-x 1 - yi), h x = -108/i. 

o o 

Below we calculate in detail the formulae for J2 and J3. The remaining relations 
can be obtained in a similar way. By changing the variables according to ([9]) and 
using (jHJ), (EJ) and Green's formula, we get (omitting the subscript 1) 

1 ff (-3 + x)(9 - Qx + x 2 + 3y 2 ) . aJ 

Jo = / / —, T7 -^-dx A dy 

54i J H < h (-3 + x-y)(-3 + x + y) 

1 I f_ o\..J_ , 1 / oN2 1 „ 3 - a; -^ 



— (b (x — 3)ydx H (p (x — 3) In dx 

18 V Jy 27 7 r ,; J 3-x + y 

1 r , 1 I ,_ 3 - X -V Jln n„ n_2 , 3n 



-J H $ In -d(27x - 9x 2 + x 

9 81 3-x + y 1 

ij -l£(27,-9^ + ^ ln |^ 
ijo - 1 jf (27x - 9x 2 + x 3 ) (^rfx - ^f^-dy 

l h ~ 8ii ( 27/2 ~ 9/3 + /4 + £ yd ^ x ~ 3)(27x ~ 9x2 + x3)) ) 

^Jo - ^( 27/ 2 - 9J 3 + J 4 - 162J a + 162 J 2 - 48 h + 5J 4 ) 
2 -((/i + 18)J -12J 2 ). 



Rewrite J 3 as the form 

J3 = J31 + J32, 

where 




J 3 i = / / (xy) x (x-y)(x + y) 



X 

(x + y) In 2(x — y) 

2/ 



dx A dy, 



6 

and 



Y. ZHAO, C. YANG AND X. CEN 



J, 




2(x — y) 2 (x + y) + (x + y)(x 2 + xy + y 2 ) 



-1^3 i „.3\ 



dx A ofa/ 




(xy) (x 6 + y 6 )dx A dy 



h<H ,o 



3 Jo. 



Note that J 31 is the integral Js(h) defined in the papers [221 IS], see §2.2 in 
or Lemma in [9J. Taking k = 1 into (jSJ), one gets the expression of J_i. It follows 
from Lemma, Lemma 1 and the appendix of [9J that 

1 -((-28-/i)/ + 24/ 2 - 12/,)). 



J: 



31 



1 / 4 1 

2 1+ 1^"6 



2 



6h 



The formula for ^(/i) follows from the expression of J31 and J32. 

The following lemma is crucial for our analysis. 
Lemma 4. [9] 27ie vector-function col(7*, I2, Iq) satisfies the following system 

( h -2 h + 6 \ 



□ 



(10) 



h 



3 -(h-6) ~(h + 9) 
{ -3 |(/i + 6) J 



Lemma 5. I' (h) < for h e (-4,0). 



Proof. Let X\(h) and X2(h) be the intersection point of and x-axis with x\(h) < 
X2(h). Note that is located in the right half-plane and has the positive (coun- 
terclockwise) orientation. It follows from (jSj) and @ that 

-1 



X 



xi(h) 



-dx 



X 



x 2 (h) 2y/hfx + (x — 3) 



-.dx + 



X2 (h) 



X 



1 ( h ) 2(-^h/x+(x-3) 2 ) 



dx 



X2Q1) 



x~ 



Define 



xi(h) y/h/x + (x - 3) 



w(h) 



-.dx < 0. 



□ 



I'M 



Let J(h) be the principal part of the displacement function 5(e, h), i.e. 
(12) J(h) = HxJ\ + 112J2 + yU3^3 + A*4^4- 
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The function J(h) can be expressed in a more convenient form by using Picar- 
Fuchs system (fTUlh Therefore, we reformulate Theorem [T] as follows. 

Proposition 6. For small e, the number of limit cycles in (j2J) that emerge from 
the period annulus is equal to the number of zeros in the interval (—4, 0) of the 
function J(h), where 



(13) 

with 
(14) 

Moreover, 

(15) 

where 



J{h) 



2(3A - 12a - IO7 + 72k) + (A - 4a + 2 7 - 8k)/i 



+ 



5184 

-A + 4(7 - 2 7 - 24k A + 67 



2592 



1296 



-1296/i 2 - 648/i 3 , 
216/i 2 - 324/13, 



a 



J'{h) 



-144/n - 468/12 - 270/13 
54/i 2 - 8I/13 + /i 4 . 

/(WO 



6/1 



4, 



7776/i(4 + /i) 



(16) /(/i) = -16A - I60-/1 + (A - 4a + 2 7 - 8k)/i 2 - 32( 7 + k/i)w. 
The corresponding integrals and function are given in Lemma and ( flTT) . 



Proof. We will prove in Corollary [23] (see Section 4 below) that J(h) = if and 
only if /i! = /i 2 = /13 = /14 = 0. Therefore, it follows from 03]) that 5(e, ft) = if 
and only if fii = /i 2 = /13 = /14 = 0. This implies that the number of limit cycles 
in (|2]) that emerge from the period annulus is equal to the number of zeros in the 
interval (—4,0) of the function J(h), defined in (|T2]) . where Jk(h), k = 1,2,3,4, 
are given in Lemma [3] 

The expression ( JTBl follows from ( IT2]) and ( ITOl . To calculate J'(h), we firstly 
derive several formulae. Differentiating both sides of system ( llOp . we have 



( h -2 h + 6 \ 
?(fc-6) ^(ft + 9) 
-3 



V 

which implies 

(17) 3/i(/i + 4) 

and 



5(^ + 6) 



( -1 \ 



J 2 



1 

4 

-- 



3 
2 
1 

2 / 



J 2 



7"" 



ft + 6 








) 
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Differentiating both sides of system (TIB"]) and then taking i fPTj) . ( TTBl into it, one 
obtains (JTSD - □ 

Remark 7. The changes ( |T4|) transform the parameters ^1,^2, ^3, ^4 into A, cr, 7, k. 
In i/ie resi of this paper, we also take A, a, 7, n as the independent parameters, 
instead of fii, i = 1,2,3,4. The function f(h) is expressed in a more convenient 
form by the new parameters. 

Without loss of generality suppose k = 1 if k 7^ 0. That is to say, k e {0, 1}. 

In the end of this section we give the asymptotic expansions of the Abelian 
integrals I^h), I 2 {h) and Io(h) near the endpoints of (—4,0). 

Lemma 8. [9j The following expansions hold near h = —4, with Iq(—A) < 0: 

h = ,. ( _4) ah±j? + w£^ + m±£t + .. 

ov ' \ 12 1296 93312 

2 oV ; V 9 1944 419904 

oV ; V 36 1944 104976 

Lemma 9. T/ie following expansions hold near h = 0: 

h = -6-lh\n 2 \h\ + --- , 
6 

27 3 

h = -y + -/iln|/i| + 3(c+l)/i + -- - , 

J = —9 + -/iln \h\ + ch-\ , 

where c is a real constant. 

Proof. The first two coefficients of every integral have been obtained in It is 
well known (see [H], or Appendix of [20J) that the asymptotic expansion of I 2 
and Iq near h = have the form <fo(h) + fi(h) In \h\, where <po(h) and (fii(h) are 
analytic at h = 0. Let 

27 3 1 
i~2 = — — + ~h\n \h\ + c 2 h H , J = ~ 9 + ~h\\i \h\ + ch H . 



Taking the above expansions into the last two equations of system f|T0|) . we get 
c 2 = 3(c+l). □ 

Using the above two lemmas, we get the following expansions for w(h). 

Corollary 10. The following expansion holds near h = —A: 

(1 9 ) w(h) = 1 + -(/* + 4) + -L(fc + 4) 2 + ^ (/1 + 4) 3 + ■ ■ • . 
The following expansion holds near h = 0: 

(20) w(h) = S + JL + .... 
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3. The geometric properties of the ratio w(h) 
In this section we study the monotonicity, convexity, etc. for the ratio w(h). 
Lemma 11. Then the ratio w(h), defined in (TlTI) . satisfies the following equation 

(21) 3h{h + 4)w' = -2w 2 + 2(h + 6)w - 2(2h + 9). 

Proof. The equation follows from system (fTTj) . □ 

Proposition 12. u>(-4) = 1, w(0) = 3, and w'{h) > for h G (-4, 0). 

Proof. The equalities w(— 4) = 1 and w(0) = 3 are obtained from (fT9|) and 
(|20|) . Since the right hand of equation (121]) is a quadratic polynomial in w and 
(2{h + 6)) 2 - 4(-2)(-2(2/i + 9)) = 4/i(/i + 4) < for h G (-4,0), we have 
-2w 2 + 2(h + 6)w - 2(2h + 9) < in the interval (-4, 0). It follows from §ZQ 
that w'(h) > for h G (-4, 0). □ 

Assume that <p(xi, x 2 ) and <p(xi, x 2 ) are two polynomials in x 1 ,x 2 - Eliminating 
the variable x« from the equations <p(xi, x 2 ) = (j>{x\, x 2 ) = 0, we get the resultant 
of ip and (ft, denoted by Resultant (<£>, 4>, x«). 

Proposition 13. For h G (-A,0),w"(h) > 0, w"'(h) > . 

Proof. Rewrite (1211) as the form 

(22) W '(h) = -^ + 2(^,-2(2^9) 
v ; v ; 3/z(/i + 4) 

By abuse of notation we denote u>' the right hand of ( 1221) . which is a function of 
two variables h and w. Then 

dw' dw' . 2(6 + h-2w)(-2h-6w + hw + 2w 2 ) 

(23) w (h) = —— + — — w - 



dh dw 9h 2 (A + h) 2 

Similarly let w" be the function of the right hand of (1231) . Using (T2"2"|) and 
yields 

(24) w'"(h) = ^ + ?fw>- AC{h > w) 



dh dw 27/i 3 (4 + /i) 3 ' 
where 

C(h,w) = -32A - 108/i -37/i 2 -4/i 3 + 2(108 - A2h + Ah 2 + tf)w 
+ (-144 + 76/i + h 2 )w 2 - 12(h - 6)w 3 - 12w 4 . 

Let 

Q = {(h,w)\-4<h<0, l<w<3}. 

It follows from Proposition [121 that (h,w(h)) G for /i G [—4,0]. 

The maximum and minimum for ((h,w) in B necessarily occur either on the 
boundary <96 of the set 9, or the points inside 9 whose coordinates satisfy 
equations 

(25) a = | = o,a = ^ = o. 
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Note 

((-4, w) = -120 - 1) 2 (19 - 8w + w 2 ), C(0, w) = -12(w - 3) 2 (w 2 + 3), 

C(h, 1) = -2(4 + h) 2 {6 + h), ((h, 3) = 2h 2 {h - 2). 

This implies that ((h, w) < in «96\{(-4, 1), (0, 3)} and C(-4, 1) = C(0, 3) = 0. 

On the other hand, we will show that ((h,w) has a unique negative minimum 
value inside G. To this end, we eliminate the variable h or w from the equations 
d2SD to get 

Resultant {( h , ( w , h) = -1024(w - 3) 2 {w -2)(w- l)\i(w), 
Resultant(Cfe,C»,^) = -QUAh 2 {h + 2)(h + A) 2 X 2{h). 

where 

XiO) = -21975 + 14660w - 1841w 2 - 912w 3 + lUw 4 , 

X2 (h) = -170856 - 4036/i - 401/i 2 + 304/r 3 + 38/i 4 . 
It is easy to obtain 

X[{w) = 2(w-2)(228w 2 -912w-3665), xW = 2(/i + 2)(76/i 2 + 304/i - 1009). 

It is readily verified that Xii w ) (resp. X2Q 1 )) nas a unique zero w = 2 (resp. 
h = —2) in the interval (1,3) (resp. (—4,0)). Thus, Xi( w ) (resp. X2{h)) has a 
maximum or minimum value at w = 2 (resp. h = —2). By direct computation we 
have xi(l) < 0, xi(2) < 0, Xl (3) < (resp. X2 (-4) < 0, X2 (-2) < 0, X2 (0) < 
0), which implies Xi( w ) < (resp. X2{h) < 0) for w G (1,3) (resp. h G (—4,0). 
Therefore, the equations f )25|) has a unique solution (—2, 2) in Q\dQ. This yields 
that ((h, w) has a minimum value ((— 2, 2) = —16. Consequently Q{h, w) < for 

(A,w)e9. 

Finally we prove ((h,w) < for (h,w) G 6\<90. We have known ((h,w) < 
for G 6. If G G\dQ and CO*,™*) = 0, then ((h,w) has a 

maximum value at (h*,w*) satisfying ( h (h*,w*) = ( w (h*,w*) = 0. By the above 
discussion we know (h*,w*) = (—2,2), which yields a contradiction. 

It follows from ((h,w) < and flMp that > for h G (-4,0). Hence 

w"{h) is an increasing function in (—4,0). By (fT9l one gets «/'(— 4) = 1/54 > 0. 
Consequently w"{h) > for h G (—4,0). This completes the proof. □ 

Corollary 14. Let 

(26) Z 1 (/i) = l + i(/i + 4), Z 2 (>0 = 3+£. 

6 2 

T/ien m t/je hw-plane, (h,w(h)) G X? ; where 

(27) D= -4 < h < 0, < w < k{h)}. 

Proof. Proposition [T2] and Proposition [TBI show that w(h) is increasing and strictly 
convex in [—4, 0). The points Si(— 4, 1) and S^O, 3) are the endpoints of the curve 
w(h) in /iw-plane. Note that u> = l\{h) is the tangents to w{h) at Si and u> = £ 2 (/i) 
is the straight line passing through both S\ and S%. The assertion follows. □ 
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4. Proof of Theorem [2] 

In this section we prove Theorem [2j 

Proposition 15. If k — 0, then J(h) has at most three zeros in (—4,0) 

Proof. If k = 7 = 0, then f{h) = (h + 4)(-4A + (A - Aa)h), which implies that 
f(h) has at most one zero in (—4, 0). It follows from (fT5l) that J'(h) has at most 
one zero. Hence, J(h) has at most two zeros in (—4,0). 

If K = 0, 7^0, then /(/i) = -327w-16A-16a/i+(A-4a + 27)/i 2 . Therefore, 
it follows from Proposition [H that f"'(h) = -32jw"'(h) ^ 0. Noting /(-4) = 0, 
this yields that f(h) has at most two zeros. Consequently, J(h) has at most three 
zeros. □ 

Now we consider the 7^ 0. To simplify the proofs, we suppose k = 1, as 

mentioned in Remark [71 By ()16p . we get 

(28) f"'(h) = -96w"'(h)p(h), 
where 

(29) p(ft) .^) + l (A + 7) . 
Lemma 16. For h G (-4,0), p'(h) < 0. 

Proof. By abuse of notation let p be the right hand of equation ( I29p . Then it 
follows from (ESI) (E2D, (E3D and (El) that 



(30) P\h) = % + ^ W '- ^Mh,w)MhM 



dh dw 2187/i 6 (4 + /i) 6 (w"'(/i)) 2 ' 
where 

Vi{h,w) = 2(/i 3 + 13h 2 + 108/i + 324) - (/i 3 - 2/i 2 + 48/i + 432) w 

-2{h 2 + 4/i - 36)w 2 , 
yjt 2 (h, w) = -2(5/i 3 + 47h 2 - 324) + (5/i 3 + 62/i 2 - 228/i - 1080)w 

-8(/i 2 - 23h - 63)w 2 - 36{h + 2)w 3 . 

In the rest of the proof let h(h), (2(h) be the linear functions given in ( 126]) and 
let V be the set defined in ( I2"7j) . Recall that Corollary IT41 shows (h,w(h)) G V. 

Firstly consider the algebraic curve *&i(h, w) = 0. Note that *&i(h, w) is a 
quadratic polynomial in w. Since —2{h 2 + Ah — 36) > 0, h(h)) = — 2(— 9 + 

/i)(4 + /i) 3 /9 > and ^(Ti, / 2 (/i)) = -h 2 (A + h) 2 < for /i G (-4, 0), there is a 
unique branch of the algebraic curve "$?i(h, w) — for (/i, to) G V\dV, defined by 

_ /r 3 - 2/i 2 + 48/i + 432 + h(h + 4)^ + 4/^ + 324 
W( } ~ 2 (-2/i 2 - 8/1 + 72) ' 
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This shows that w(—A) = 1, w'(—4) = 1/6, w"(— 4) = 13/162. By direct com- 
putation, one gets for w G (1, 3), 

Resultant f * la ^-i, /» ] = -746496(w - 3) 2 (w - 1) 2 (69 - 20w + 5w 2 ) < 0, 
\ ow J 

which implies that ^i,,, ^ on the algebraic curve *&i(h, w) = 0. 

Now we are going to show that the curves w(h) and w(h) do not intersect in 
V\dV. Let 

v(h) = w(h) — w(h). 

By direct computation, v(h) = —5(h + 4) 2 /162 + • • • < as h — > —4, and 
v(h) = 6/ In \h\ + • • • < as /i — )■ 0. Hence, v(h) have at least two zeros in 
(—4,0), counted with their multiplicities, if the curves w{h) and w{h) intersect 
in T>\dT>. Assume that hi and h 2 are two zeros of v(h) with hi < h 2 if possible. 
This yields that v'(h) has at least one zero h* such that h* G (hi, h 2 ) if hi ^ h 2 , 
or h* = hi = h 2 . That is, the following equations have at least one solution for 
h G (-4,0): 

31 w) = 0, w (h) + — = , , ,x. Tf = 0, 

Wi w 3/i(/i + 4)Wi^ 

where it/(/i) is defined in fT22|) . and 

= 2(ll/i 4 + 119/i 3 + 714/i 2 + 2592/i + 3888) - (ll/i 4 + 16/r 5 + 32/i 2 

+2304/i + 7776)w - 6(3/i 3 + 26/i 2 - 16h - 432)w 2 

+8(/i 2 + 4/i-36)w 3 . 

However, 

Resultant^, *i, w) = -466560/i 5 (/i + 4) 5 (h 2 + Ah - 36) < 0, fte (-4, 0). 

Consequently, the equations (l3Tj) has no solution for /i G (—4,0), which yields 
contradiction. Thus, the curve w(h) does not intersect w(h) in the set V\dV and 
w(h) < w(h) for /i G (—4, 0), which means *f?i(h, w(h)) ^ 0. Since li(h) < w(h) < 
w(h) and ^i(h,h(h)) = -2(-9 + /i)(4 + /i) 3 /9 > 0, one obtain ^(/i, > 
for h G (-4,0). 

Secondly we consider the algebraic curve ^f 2 (h,w) = 0. Note that ty 2 (h,w) is 
a cubic polynomial of w for /t 7^ —2. If h ^ —2, then ^(Z?-, w) = has at most 
three branches. Since for h G (—4,0), 

* 2 (h, 1) = " 5/ ^ + 4) 2 > 0, # 2 (/i, ^(/i)) = - 3)(h + 4) 3 < 0, 

9 

h(h)) = -4h 2 (h + 4) 2 < 0, V 2 (h, 3) = 5h 2 (h + 4) > 0, 

we conclude that the algebraic curve ty 2 (h,w) = has no branch in T>\&D, one 
branch in 1 < w < li(h) and one branch in l 2 (h) < w < 3 respectively, provided 
that h G (-4,0) and h ^ -2. Therefore, $ 2 (h,w) < for (h,w) G V\dV, 
h^-2. 
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If h = —2, then ^ 2 ( — 2, = 8(44 — 52w + 13w 2 ). Direct computation shows 
that \&2( — 2,w) = has two zeros at 

wl = A ^13 _ y/^fj ~ 1.21554, = ^ (l3 + v 7 ^) « 2.78446. 

On the other hand, the straight line h = —2 intersects the boundary of the 
set V at two points (-2,4/3) and (-2,2), which implies that 4/3 < w* < 2 if 
(-2,w*) G V. This gives (-2,<) g V, i = 1,2. Hence, * 2 (-2,tu) < for 
(-2,w) G 2?\9P. 

Summing the above discussions, we get < from ( |30l) . □ 

Proposition 17. Ifn = l, 7 G (—00,— 26/7] U [0,+oo) 7 t/jen J(/i) /ias at most 
three zeros in (—4,0). 

Proof. It follows from (JTSJ and flU} that 

(32) p(-4) = 1(7 7 + 26), p(0) = lim p(/i) = i 7 . 

Zl h-*0 d 

It follows from Lemma [16] that p(/i) > if 7 G [0, +00) and p(h) < if 7 G 
(-00, -26/7], G (-4,0), which implies that f"(h) 7^ for h G (-4,0) by Q2S} 
and Proposition [131 Hence, f ,f (h) has at most one zero. Noting /(— 4) = 0, /(7i) 
has at most two zeros in (—4,0). One obtains by using f [T5l) that J(/i) has at 
most three zeros. □ 

The following lemma gives the asymptotic properties for f{h) and J(h). 

Lemma 18. Let k = 1. T/ie following identities hold: 

/M) = 0, /'(-4) = -|(3A - 6,7 + 87 - 20), 

J(-4) = ^(-4), J-M) = -— /< (-4«(-4). 
T/ie following expansions hold near h = 0: 

J^L\n 2 \h\ + --- , Jv 

124416 11 ' J v ' J w ln|/i| 



(34) J(h) = J^- In 2 !/»! + ..., f(h) = /(0) - + 



/(0) = -16(A + 6 7 ). 

Proof. The results follows from (1131) . ( II 6p . (j!9p . ( I20p . Lemma [8] and Lemma M D 

Proposition 19. Suppose « = 1, 7 G (-26/7,0) and h G (-4,0). 1/ /(0) > 0, 
i/ien J(/i) /ias at most tnree zeros; If /(0) < ; inen J(/i) aas at most two zeros. 

Proof. If /c = 1, 7 G (—26/7,0), then it follows from (|32|) and Lemma [TBI that 
p(/t) has a unique zero in (—4, 0), which implies that f"(h) has at most two zeros 
by (!28|) . Hence f'(h) has at most three zeros. Noting /(— 4) = 0, has at 
most three zeros in (—4,0). This yields that J'(h) has at most three zeros by 
(|T5|) . Therefore J(h) has at most four zeros in (—4,0). 

If 7 < 0, then it follows from (1331) and Lemma [5] that J(— 4) < 0. 
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We split the proof into three cases. 
Case 1. /(0) > 0. 

It follows from (134")) that J(h) > as h — > 0. The number of zeros of J(h) is 
odd. Hence, J(h) has at most three zero in (—4,0). 
Case 2. f(0) = 0. 

Since f"(h) has at most two zeros and /(— 4) = /(0) = 0, f(h) has at most 
two zeros in (—4,0). It follows from (fT5|) that J'(h) has at most two zeros. This 
implies that J(h) has at most three zeros in (—4,0). 

Case 3. f(0) < 0. 

By (J33D J(h) < as h ->■ 0. Noting J(— 4) < 0, the number of zeros of J{h) is 
even. 

If /'(— 4) < 0, then the number of zeros of f(h) is even. Therefore f(h) has at 
most two zeros in (—4,0). This means that J(h) has at most three zeros. Since 
the number of zeros of J(h) is even, J(h) has at most two zeros. 

If /'(— 4) = 0, then f'(h) has at most two zeros. Since /(— 4) = 0, f(h) has 
at most two zeros. Consequently, J(h) has at most three zeros. Noting that the 
number of zeros of J(h) is even, J(h) has at most two zeros. 

It follows from Lemma M that f'(h) = 1927/(Mn 2 \h\) + ■ ■ ■ > as h -» 0. If 
/'(— 4) > 0, then the number of zeros of f'(h) is even. Since f'(h) has at most 
three zeros, f'(h) has at most two zeros in (—4, 0), which shows that f(h) has at 
most two zeros in (—4,0). Therefore it follows from (Tl5l) that J'(h) has at most 
two zeros, and hence J(h) has at most three zeros. Since we have shown that 
the number of zeros of J(h) is an even number, J(h) has at most two zeros in 
(-4,0). 

The proof is finished. □ 

By the above discussions we know that J(h) has at most three zeros in (—4, 0). 
In what follows we are going to show that there exist the parameters //j,z = 
1, 2, 3, 4, such that J(h) has exactly three zeros in (—4, 0). To do this we introduce 
the following definitions (see for instance [6l [TT]). 

Let /o, /i, • • • , fn-i be analytic functions on an open interval L of R. (/ , /i, • • • , f n -i) 
is a Chebyshev system on L if any nontrivial linear combination 

Ao/o(a?) + AiA(x) H h A n _i/ n _i(x) 

has at most n — 1 isolated zeros on L. 

{fo, fi, ■ ■ ■ , fn-i) is a complete Chebyshev system on L if (/ , /i, • • • , fk-i) is a 
Chebyshev system for all k = 1, 2, • • • , n. 

(/o) A) - '' j /n-i) is an extended complete Chebyshev system (in short, ECT- 
system) on L if, for all = 1, 2, • • • , n, any nontrivial linear combination 

(35) Ao/o(aO + Ai/i(ar) H h \ k -ifk-i{x) 

has at most fc — 1 isolated zeros on L counted with multiplicities. 

Remark 20. // {fo,fi,--- , fn-i) is an ECT-system on L, then for each k = 
1, ...,n, then there exists a linear combination (|35|) with exactly k — 1 simple zeros 
on L (see for instance Remark 3.7 in [7]j. 
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Lemma 21. [6J [TT](/o, fi,"' > fn-i) an ECT-system on L if, and only if, for 
each k = 1, 2, • • ■ , n, 



./«» 
./,'• 


(s) 
(*) 


• 


■■ /*-i(ar) 

•• fUV) 


.(&- 




15 (*) 







?0. 



Proposition 22. There is a real number 6 G (—4,0] such that (Ji, J2, J3, J4) zs 
an ECT-system in (—4,6). 

Proof. In this proof we use the notation introduced in Lemma I2T1 Let /fe(/t) = 
Jk+i, k = 0,1,2,3. By Lemma [3] and Lemma we have Ai = J\{h) < 0. It 
follows from Lemma [8] that the following expansions hold near h = —4, with 
a Q = 7'(-4) < 0: 

A * - a »(-i^ (ft + 4)2 + "')' 



A > = a »(5S784 (,! + 4)3 + " 

a 4 = osC ! 47 " i+4) + . 

V 6377292 229582512 

Therefore there is a real number b G (—4, 0] such that Aj(/i) < 0, 2 = 1, 2, 3, 4, for 
/i G (—4,6). The assertion follows from Lemma [2TT □ 

Corollary 23. J(ti) = if and only if fa — 0, i = 1, 2, 3, 4. 

Proof. Suppose that J(h) have the form as f|T2l . If J{h) = 0, then J'{h) = 
J"{h) = J"'(h) = 0. In the proof of Proposition (122]) . we have shown that A4 < 
as h -> -4, which implies that J(/i) = J'(/i) = J"(/i) = J'"(/i) = if and only if 
/ij = 0, i — 1, 2, 3, 4. This completes the proof. □ 

Now we can prove the main results of this paper. 

Proof of Theorem^ By Remark [?l Proposition [T5l Proposition [TTl and Proposi- 
tion [191 J(h) has at most three zeros in (—4,0). Proposition [221 and Remark [201 
show that there exist fa,i = 1,2,3,4, such that J(h) has exactly three zeros in 
(—4,6), 6 G (—4,0]. The theorem follows from Proposition [6j □ 
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